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We derive the explicit values of all regularization parameters (RP) for a scalar particle in an 
arbitrary geodesic orbit around a Schwarzschild black hole. These RP are required within the 
previously introduced mode-sum method, for calculating the local self-force acting on the particle. 
In this method one first calculates the (finite) contribution to the self-force due to each individual 
multipole mode of the particle's field, and then applies a certain regularization procedure to the 
' mode sum, involving the RP. The explicit values of the RP were presented in a recent Letter [Phys. 

■ Rev. Lett. 88, 091101 (2002)]. Here we give full details of the RP derivation in the scalar case. 
' The calculation of the RP in the electromagnetic and gravitational cases will be discussed in an 
, accompanying paper. 

o ■ 

; I. INTRODUCTION 

^ . The space-based gravitational wave detector LISA (Laser Interferometer Space Antenna), scheduled for launch 

■ around 2011 Hi, will open up a window for the low frequency band below IHz, allowing access to a variety of black 
0^ . hole sources. As one of its main targets, LISA is expected to detect the outburst of gravitational radiation emitted 

■ during the capture of a compact star by a supermassive black hole — a 10^~^ solar masses black hole of the kind now 
believed to reside in the cores of many galaxies, including our own [||. Designing accurate gravitational waveform 

' templates for this type of astrophysical events requires an accurate knowledge of the orbital evolution, including 
the effect of radiation reaction. The evolution of such extreme mass-ratio systems can be modeled by considering 
a pointlike test particle moving in the fixed gravitational field of a black hole. One then addresses the question of 
I the local self force acting on this particle. (In special cases, one may study the orbital evolution under radiation 
I , reaction using global energy- momentum balance techniques [||. Such techniques, however, appear insufficient when 
dealing with the astrophysically realistic case of non-equatorial eccentric orbits in Kerr spacetime.) A general formal 
scheme, well physically motivated, for calculating the gravitational self-force in curve spacetime was worked out by 
Mino, Sasaki, and Tanaka (MST) Q and, independently, by Quinn and Wald (QW) Recently, the two groups of 
Barack and Ori (BO) and Mino, Nakano, and Sasaki (MNS), have reported on a practical method for implementing 
this formal scheme |^, allowing, for the first time, actual calculations of the gravitational self- force for any geodesic 
d ' orbit around a Schwarzschild black hole. The purpose of the present paper (together with the one accompanying it 
[0 ) is to provide a full account of the method and results reported in ||] . 

The notion of self-forces is briefly described as follows. Consider a pointlike particle carrying a charge q, which may 
represent here a scalar charge, an electric charge, or a mass. The particle is assumed to move freely in the curved 
background of a black hole with mass M ^ q. In the limit q —> 0, such a particle is known to move along a geodesic 
of the background geometry. However, when endowed with a finite charge (or mass), the particle no longer traces a 
background geodesic, as a result of interaction with its own field. The finite-charge (or finite-mass) correction to the 
particle's motion is then described in terms of a "self-force" : Treating the particle's field as a linear perturbation on 
the fixed black hole background, the particle's equation of motion is written as 



o 



- 

X 



/^aa=Fr, (1) 

where /i is the particle's mass, Oq denotes its (covariant) four-acceleration, and oc 0{q^) describes the leading- 
order self-force effect. [In the gravitational case, the four-acceleration, as well as the ("fictitious") self-force, are 
defined through a mapping of the particle's worldline into a trajectory in the background spacetime — see Ref. |^.] 
The formal construction of was described by MST Q| and QW in the gravitational case, by DeWitt and 
Brehme and QW |^ in the electromagnetic case, and by Quinn in the scalar case. In all cases, the self-force 
is constructed through 
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F^"'* = Km F^'^'^x) + trivial local terms, (2) 

X — >z 

where z represents a point on the particle's worldline where the self-force is being evaluated, x is a point in the 
neighborhood of z, and the local terms are given explicitly in (they include the Abraham-Lorentz-Dirac force 

in the scalar and electromegnetic cases). The quantity F^^^^{x), the "tail" force, is a nonlocal contribution to the 
self-force, whose occurrence reflects the essential nonlocal nature of the radiation reaction effect in curved spacetime: 
waves emitted by the particle may backscatter off spacetime curvature and later interact back with their emitter. The 
tail force is formally constructed through a worldline integral as P,|5|,|l0| 

F^-i(.x) = hm q r 'w^G[x,z(T)]dT, (3) 

where t is the proper time along the particle's worldline, tq is the value of r at the intersection of the worldline 
with the past light cone of x, G symbolizes a Green's function for the particle's field, and Vq is a certain first-order 
differential operator acting on G [the explicit form of V^, as well as the type of the Green's function — whether a 
bi-scalar, a bi- vector or a bi-tensor — depend on the case considered — whether scalar, electromagnetic, or gravitational 
(see [pp|JlO[| for details)]. Notably, in the gravitational case, the tail force constitutes the sole contribution to the 
self- force (as long as geodesies in vacuum background are considered) . It is the actual evaluation of the tail part that 
has rendered practical calculations of the self-force most challenging. 
It is instructive (and later useful) to write Eq. (^ in the form 

Fi^^{x)^F'^\x)-Ff\xl (4) 

where F^^^{x) and F^"{x), the "full" and "direct" forces, are the quantities constructed by replacing the integral in 
Eq. d) with /^'^' and ^^1, respectively. The "full" force F^"" is directly obtained from the particle's "full" field 

by acting with gVa [for the scalar case, e.g., see Eq. ( |l5| ) below]. The "direct" force F^" is the "divergent piece" to 
be removed, which is associated with the instantaneous effect of waves propagating directly along the particle's light 
cone. Note that the "tail" force is hence attributed to waves scattered inside the particle's past light-cone. 

A direct implementation of the MST and QW scheme for calculating the self-force in weak field was introduced 
recently by Pfenning and Poisson [ pT| . To allow a practical implementation of this formal scheme for strong-field or- 
bits, BO devised a multipole-mode decomposition method, relying directly on MST and QW's formal result (gh. BO's 
mode-sum method was formulated first for the scalar self- force [ p^ , and later for the gravitational self- force [ p^ . This 
method has been tested and fully implemented for calculating the scalar self force in several cases [ p^|jl5[| . The mode 
sum scheme (which we review in the next section) is based on decomposing the tail force into individual multipole- 
mode contributions, relating these contributions to the "full force" modes — which are accessible to standard numerical 
analysis — and then summing over the mode contributions, subject to a certain regularization procedure. This proce- 
dure requires knowledge of certain analytic parameters, the "regularization parameters" (RP), whose values depend 
on the orbit under consideration. The RP values were derived previously for a few special orbits in Schwarzschild 
spacetime: for radial and circular orbits in the scalar case [ p^ , and for radial trajectories in the gravitational case 
[ p^ , p^ . These (rather cumbersome) calculations were carried out through a special local perturbative expansion of 
the Green's function's multipole modes, relying directly on the integral formula (||). 

In this paper we present a different approach for the calculation of the RP, based on a direct multipole decomposition 
of the "direct" piece of force. This new approach (already outlined in [^) allowed a rather convenient calculation 
of all RP values for a general geodesic orbit in Schwarzschild spacetime, as we describe in this paper. In particular, 
it provided an independent verification for the RP values in the special cases considered previously (using the l- 
mode Green's function analysis as mentioned above). Two variants of the new calculation method were worked out 
independently by the two groups of BO and MNS, yielding the same RP values The calculation by MNS has been 
reported in [0. This paper presents full details of the RP derivation by BO. 

In its basis, the calculation method presented here is applicable to all three sorts of self-forces: scalar, electromag- 
netic, and gravitational. We find it most instructive to concentrate first on the scalar toy model. This 
model captures the essential parts of the calculation technique, while avoiding several complexities and delicate issues 
that show up in the gravitational and electromagnetic cases. In this paper we thus focus on the scalar model, leaving 
the treatment of the gravitational and electromagnetic cases to an accompanying paper. 

It should be commented that other approaches for the calculation of the self-force, not directly relying on the MST 
and QW formal scheme, were also suggested recently. Lousto introduced an approach also based on a multipole 
decomposition but employing a proposed zeta-function regularization scheme. Other methods were proposed by 
Nakano and Sasaki [|l^, Detweiler [M, and Detweiler and Whiting [|2l|]. 
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The paper is arranged as follows: In Sec. || we review the mode-sum method, and define the regularization 
parameters. The scalar toy model to be considered in this paper is introduced in Sec. III. The expression for 
the "direct" part of the scalar force is introduced and processed in Sec. IV, and is being formally decomposed into 
modes in Sec. 0. We then prepare for the calculation of t he R P by introducing a useful coordinate system, in Sec 



VII 



wher 



e, 



through an investigation of the d irect force's 
Sec. 



0. The main part of our calculation is contained in Sec 

multipolc modes, we obtain all RP values for a general trajectory in Schwarzschild spacetime 
the RP values, and Sec. [X provides some concluding remarks. 

Throughout this paper we use geometrized units (with G = c = 1), and metric signature — 
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summarizes 



II. REVIEWING THE MODE-SUM APPROACH 

The mode-sum method was introduced in Ref. jl^ for the scalar self force, and in Ref. [|l3| for the gravitational 
self- force. Here we review it using a slightly different perspective (and notation). 

In the mode-sum scheme, one first formally expands all three quantities F^'^^^{x), F^"''(x) and F^"{x) appearing in 
Eq. (^ into multipole Z-modes as 

oc oo oo 

1=0 1=0 1=0 

(where, recall, x represents an off-worldline point in the neighborhood of the self- force evaluation point z). Here, 

are the quantities obtained by summing over all azimuthal numbers m (and, in the 
gravitational case, also over all ten tensor harmonics), for a given multipole number /. An important benefit of the 
multipole decomposition is the fact that, whereas i^^"" and F^" both diverge a,t x ^ z, their individual modes attain 
finite values even at the particle's location (though they are usually found to be discontinuous there). Applying the 
multipole decomposition to Eq. (Q) we obtain 

^^taii);(^) = Fif"")'(:E) - Fi^^'^^^ix). (6) 

Considering now MST and QW's expression for the self-force, Eq. (H), we have 

pscU ^ ^_tail(^ = ^) = E = ^) (7) 

(we hereafter ignore the trivial local terms and focus on the tail contribution) . Note that since the tail force F^^"(j;) 
is regular at the particle's location z QJ^, one gets F^'^'^ by just evaluating the tail force a.t x — z. We can then write, 
using Eq. (|), 

^ ^ X — >Z ^ ^ \_X — ''Z X — J 

I I 

where the direction of the limit a; — > z is considered as prescribed. It is important to note here that each of the two 
limits YiTHx^z Fa"^^^^^ {x) and \\Tax=,z F^''^^^ {x) is, in general, directional dependent. This, however, does not pose a 
problem (and the third equality in the above chain of equalities is valid) if the direction of the limit is prescribed: 
one then only has to make sure that the two limits of the full and direct forces are taken in a consistent manner (i.e., 
from the same direction). 

In the last expression of Eq. (^), the sum over ^-modes is guaranteed to converge [as F^^^{x) is a regular function]. 
However, the individual sums over the full-force modes and over the direct-force modes usually diverge. Suppose now 

that one could construct a function /i^ that would make the sum lim^.^^ ^"^^"''''(2:) — h}^ convergent. Then, we 

would have [continuing the chain of equalities (0)] 



I 



P'^'' = E [(li- - h'o) - (lim f^r^)'(x) - hi] 

'hm F(*"")'(x) - /iL) - E (1™ - ^0 ■ (9) 

^X — ''Z / ^ — ^ \x — >z / 



I 



In principle, the "regularization function" is to be obtained by exploring the behaviour of the full-force modes at 
large I. However, this function can also be deduced by analyzing the large I behavior of the local quantity F^"^^ — a 
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task accessible to analytic treatment. In all cases considered so far, the function /i^ was found to have the general 
form 

hi = A^L + B„ + Co^/L, (10) 

with L = I + 1/2, and where A^, 3^, and Cq, are Z-independent coefficients whose values depend on the details of the 
trajectory under consideration. Defining 



OG 

D^ = y] [lim Fi''"^'{x) - A^L - - C^/l] , (11) 



/=0 



we finally get from Eq. 



K"'' = E - AaL - B„ - C^/l] - D„. (12) 

Eq. ( |l^ ) constitutes the basic formula for constructing the self-force through the mode-sum method. The four 
quantities Aa, Bey, Ca, and Da are called the "regularization parameters" (RP). The full modes F^^^, recall, are 
directly obtained from the "full" field modes [see Eq. (|l^) below for the construction of the full force in the scalar 
case], which, in turn, are calculated using standard numerical techniques. Eq. (|l^) thus describes a practical scheme 
for constructing the self force, given the values of the RP. 

In this paper (dealing with the scalar self-force) and in the accompanying paper (dealing with the gravitational and 
electromagnetic self-forces) we derive the values of all RP needed for implementing Eq. (^2|) for any geodesic orbit in 
Schwarzschild spacetime. 



III. SCALAR TOY-MODEL 



We consider a particle of a scalar charge q, moving freely in the vacuum exterior of a Schwarzschild black hole with 
mass M q. In the lack of self- force, the particle moves along a geodesic z'^{t) with specific energy and angular 
momentum parameters £ and £, respectively. We shall consider the self-force acting on the particle at a point along 
its worldline which we denote by z = {to, tq, ^o, ^o) (where t, r, 9, (p are the standard Schwarzschild coordinates). Let 
also X = {t, r, 0, (p) denote a point in the close neighborhood of z. 

The particle induces a scalar field $^""(x), which we shall treat as a linear perturbation over the fixed Schwarzschild 
background. In our model, the field <I>^""(a;) is assumed to satisfy the (minimally-coupled) Klein-Gordon equation 

where a semicolon denotes covariant differentiation with respect to the background geometry, and the scalar charge 
density is given by 

/"OO 

p{x)=q S*[x - z{r)]{-g)-'/^dT (14) 



{g being the metric determinant). We now define the "full force" as the vector field 

Fi-'\x) EE q^%^ (15) 

Note that both the full field <I>^"''(a:) and the full force F^^^{x) obviously diverge on the worldline, but are otherwise 
well defined. 

The force definition ( [l5| ) complies with Quinn's definition |]l0||. It differs from the expression used by MNS, which 
involves a spatial projection of the scalar force [see Eq. (1.3) in ]17||]. We prefer to adopt here the force definition ( p^ ) 
for several reasons: (i) It is a simpler definition, which nevertheless serves as an effective toy-model for the realistic 
gravitational case, (ii) It avoids the need to consider an off-worldline extension of the four-velocity, as necessary for 
defining the spatially-projected force, (iii) The force model ( p^ is naturally derived from a Lagrangian formalism, 
and is hence consistent with global stress-energy conservation — unlike the spatially projected force |10| , |22|]. 

Finally, we introduce the notions of the "direct" field and the "tail" field = <I>f"" - (sce ^^|jl7| ), from 
which the direct and tail forces are derived by 

F^"ix) = q^'^t', F'a^'\^) = Q^'f- (16) 

Recall that the "direct" field is the part of the scalar field propagated directly along the particle's light cone, while 
the "tail" part is associated with reflections of the field inside the light cone. 
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IV. DIRECT FORCE: PRELIMINARIES 



The form of the dhect scalar field was worked out by MNS (see also some preliminary results in by 
studying the Hadamard expansion of the field equation. Let e{x) denote the spatial geodesic distance from the point 
X to the geodesic z{t) (i.e., the length of the short geodesic section connecting x to the worldline and normal to it), 
and let = x'^ — z''. Then, the direct scalar field obtained by MNS can be written in the form 

$dir(^) ^ i/M-,. const, (17) 

where / is a regular function of dx (and z) satisfying 

/ = 1 + 0(^2) (18) 

(the explicit form of / will not be needed in the analysis below). Introducing the squared geodesic distance S{6x) = e^, 
the direct scalar force is then given by 

F^"{x) = q'f%^ = [laS~'/^ - (/72)^-3/25„] . (19) 

Consider now the Taylor expansion of the function S{Sx) about Sx = 0. We write this expansion as 

S = Sa + Si + S2 + --- (20) 

where Sq,Si, . . . represent terms of homogeneous orders 6x'^, Sx^, . . ., respectively. Note that this decomposition of S 
is no longer covariant, and the individual terms Sn will depend on the choice of coordinate system. Below we shall 
need only the two leading terms, 5*0 and 5*1, which we obtain in Appendix |A|. We find 

So^igl, + u^u,)5x^'5x'', (21a) 

^1 = (uA^r^o + 5°/3,V2) Sx''Sx^6x\ (21b) 

where m" = dx'^/dT is the four- velocity at z, and and denote, respectively, the connection coefficients and 
metric functions evaluated at 5x = (namely, at a; = 2). Substituting Eqs. (^ and ( |2l| ) in Eq. (pj|), we now obtain 
a Taylor expansion for the direct force, which we may expresse as 



E-3p(l) , g-5p(4) , £-7p(7) 



(22) 



Here, eo = Sq^^, and Pa^^ denote terms of homogeneous order 0((5x"). Note that the term of the form cx e^^Sx in 

Eq. ( pj| ) can be written as cx e^'^{e^Sx) oc e~''Sx'^ and then be absorbed in the term e^'^Pa^. Similarly, terms of the 

form C!C e~^6x-^ may be expressed as cx e~^Sx'^ and be absorbed in Eq -^Pi'^\ and so on. Note also that the three terms 
presented in Eq. ( p^ ) are of orders Sx~-^, Sx~^ , and Sx^, respectively. The three dots (• • •) in that equation represent 



terms that vanish in the limit Sx ^ [like, e.g., Cg esc 0{Sx)]. In the following analysis we shall need the explicit 



values of only Pi^"* and Pq^\ which are given by 

Pi'^ - -^^o,a, (23a) 

Pi^'^ - -^5o5i.„ + ^So.cSi. (23b) 

Now, in constructing the self-force, one is merely concerned with the behavior of the direct force at a; — > z — see, 
e.g, Eq. (||). Thus, the terms represented by the three dots (• • •) in Eq. (|2^), which vanish in the limit x ^ z, are 
irrelevant for calculating the self-force, and may be ignored in our analysis. We hence introduce a "revised" version 
of the direct-force by omitting these terms (retaining, though, the notation F^"): 



(24) 



p{A) , p{B) , p(C) 

where 

-3p(l) p(B)^,-5p(4) p(C)^ -7p(7) (or.) 

Note that this splitting of i^^"^ holds for any choice of coordinates x^ which are sufficiently regular in the neighborhood 
of z (though the coefficients Pi""* will depend on the choice of coordinates). 
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V. MULTIPOLE DECOMPOSITION 



Next, we consider the multipole decomposition of F^". Let 

^^'''■(^)=E^^(^'*)^""(^'^)' (26) 

Ini 

where Y^'^{9, (p) are spherical harmonics. We denote by F^^ the total ^mode contribution to the direct force at z: 

,lim Vi^'"(r,to)>^""(eo,^o). (27) 

or^O^ — 
m 

Note that F^^ = lim^^^ Fi'^"^\x) [as in Eq. (||), e.g.], where the direction of the limit is exphcitly specified such 
that X approaches z "from the radial direction". The ± sign corresponds to the two possible radial limits, r 
or rather r — )■ Tq". This choice, of taking the radial limit, appears most convenient in our multipole-mode scheme. In 
particular, it is most easily implemented in the (numerical) calculation of the full force modes (recall that the limit 
X — !■ z of both the direct and full forces must be taken from the same direction) . 



Eq. (27) is invariant under rotation in the subspace of angular coordinates 9, ip. We take advantage of this property, 
and re-define the angular coordinates such that z is located at the pole, i.e., 6*0 = 0. Due to angular-momentum 
conservation, the particle is now confined to move on a plane of constant if, which we take as = 0, tt (the value of 
the (p coordinate is fixed along the particle's trajectory, apart from a "jump" at the two poles 6 = 0, tt). The particle's 
four-velocity now satisfies u''' = 0. 

The above setup is beneficial in that the Z-mode F\.^ is now composed of only the axially-symmetric m = 
harmonic: RecaU that F'™ vanishes at 61 = for any m ^ 0, and ^'■'"=0(61 = 0) = [L/{2tt)]^/'^Pi{1), where P/(cos6') 



is the Legendre polynomial and P(l) = 1. Consequently, we find from Eq. {U7\ 



= \im [L/{2n)]'/'F^"^='{r,to). (28) 



The mode F^.'»=" is given by the integral 

Fi™=0(r,t) = J F^"\r,t,9,ip)[Y'^"^-rdn 



[L/(2^)]i/2 F^"ir,t,9,^)Piicose)dn, (29) 



where dfl = d cos 6 d(p and the asterisk denotes complex conjugation. Combining Eqs. ( p8|) and ( P9|) we finally obtain 
the following integral expression for the total Z-mode direct force: 

FL= lim^ f F^'^{r,to,0,^)Pi{cose)dn. (30) 



VI. REGULAR COORDINATE SYSTEM 

The coordinate system (t, r, 9, ip) is singular a,t 9 = Oq = 0. This singularity makes the expansion (|2^), ( ]2l| ) inappli- 
cable in these coordinates. To overcome this difficulty, we introduce the two "locally-cartesian angular coordinates" 

x = p{9) cos If, y = p{9) sin ip, (31) 

where p{9) is a sufficiently regular, odd function of 9, admitting the expansion 

p{9)=9 + pi9'' + p29^ + ■■■. (32) 

For later convenience we shall also demand that p{9) grows monotonously within the entire domain < 9 < t:, such 
that p{9) is invertible. [An obvious natural choice would he p = 6; however, later we shall make the specific choice 
p{9) = 2sin(0/2) which will simplify our calculations.] 
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Using the relations dp/dO = 1 + and p^/sin^ 9 = \ + p^/i2(p^) [easily followed from the above definition 

of where hi and /12 are both regular functions of p^, one finds that the contravariant components of the metric 

tensor now take the form 

gyy ^r-^{l + y^hi+x^h2), 

g-y = r-\hi~h2)xy. (33) 

The point z is located at x = y = 0. The above tensor g"'^ is perfectly regular in the neighborhood of this point — and 
so is the covariant metric gaf3- In the particle's location itself, a; ~ y = 0, the line clement takes the simple form 

9L=9ly^rl g% = Q (34) 

[along with = -(1 - 2Af/ro) and = (1 - 2M/ro)"i]. 

Note that the particle's geodesic is confined to y = and correspondingly = 0. Also, since z is located at 
a; = 2/ = 0, we have 5x^ — x, Sx^ = y. Finally, we comment that the particle's angular momentum is given as C = 
evaluated at z (but note that is not conserved along the geodesic). 



VII. INVESTIGATING THE Z-MODE OF THE DIRECT FORCE 



A. Is the X ^ z limit interchangeable with the Legendre integral? 



We now explore in more detail the Z-mode direct force Fj-q' based on the integral formula (30). Recalling that the 



direct force itself is composed of three terms, Eq. (E4h, we write 



±a 



p{A)l , p{B)l , p{C)l 



(35) 



where Fa {W standing for A, B, or C) denotes the contribution to F\.^ [through Eq. (|30|)] due to the term Fa 
of the direct force : 



(W) 



±a 



— lim 



L 



Sr^o± 2tt 



Fi'^\r,to,9,if)Pi{cose) dn. 



(36) 



Recall that the various terms F^'' are given in Eq. (|2 



The task of evaluating the various contributions would be much simplified could we interchange the limit 

Sr — > 0^ and the integration in Eq. (^6[). Is such an interchange allowed? In Appendix ^ we address this question, 
and show that interchanging the — > 0^ limit and the Legendre integral is indeed allowed for = _B and W = C; 
however, as evident from the explicit calculation below, such an interchange is no valid for W = A. Here we present a 
heuristic argument suggesting why the interchange is valid for W = B,C, and why it might fail for W = A. A sketch 
of a mathematical proof is provided in Appendix |b[ 

For a given small separation 6x = x — z, assume that all components of Sx"' scale as Sr (we assume Sr ^ 0). Since 
eo then scales like Sr too, we find that the various terms F^^-* scale as 



p(4) 



piC) _ ,-7p(7) 



(37) 



where the proportion coefficients only depend on the "direction" of Sx'^ (i-e., on the ratios between its various 
components). To consider the interchangeability of the limit and integral in Eq. (^6|), one is mainly concerned with 
the contribution to the integral from small x,y values (i.e., from the immediate neighborhood of the integrand's 
singular point z). To find out how this small piece of integral scales with (5r, we consider the small integration area 
around z, in the xy plain, defined by p = (x^ + y^Y^"^ < Sr (for a given Sr ^ 0). Observing that this integration 
area scales like Sr^ and relying on the scale relations (|37|), one finds that this small-^r contribution to the integral 



like Sr^ for f!^^\ and like Sr^ for f!^^K Namely, upon taking the limit Sr 0, the small-5r 



scales like for f!^^'' , 

piece of integration vanishes for W = B,C, but not for W ~ A. This suggests that we may interchange the limit and 
integration for W = B,C, but not for W ^ A. A more rigorous mathematical treatment implies that this is indeed 
the case (see Appendix ^ . 
We are thus allowed to write 
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F^^^^^^ - ^1 Fi^^^\r^,t,,x,y)Pi{cose)d^. (38) 
However, for W = A we must use the original expression: 

f'^^^ = Wui [ Fi^\rM,x,y)Pi{cose)dn. (39) 

i5r— >0± ZlT J 

For later convenience we give here explicitly the form of ' for r = ro,t — to. We have 

Fi^^ = e,'Pi^^ (x, y), Fi^) = r,' P^J^ (x, y), (40) 

where P^\x^ y) is a polynomial of homogeneous order n in a; and y, and eg is the reduction of eq to 6r = St = 0: We 
find, recalling = + y^ and Mj, — 0, 

eo = (»-oP + ^xs; ) • (41) 
Note that eo is an even function of both x and y — a fact that will play a crucial role in the analysis below. 



B. Calculating fI''^' 



Let us first evaluate the term Fa ■ We observe that the integrand in Eq. (38) is composed of three factors: 



Eq ^ ^ ^« i^iV) X -Pi(cos0). Since eg and cos9{p) are even function of both x and y, then so are the factors Cq ^ 
and P;(cos6'). However, each of the eight terms of Pa'^ (x, y) (proportional to a;^y°, x^y, . . . , x°j/^) is of odd power in 



either x or y. Hence, the overall integrand in Eq. ( |38[ ) is composed only of terms which are odd in either x or y. As 
a consequence, the integral is found to vanish identically, yielding 

i^^C^) = 0. (42) 



C. Calculating 

We next turn to consider Fq*-^-* . The integrand in Eq. ( p8| ) now takes the form e^^^ x Pi^^ {x,y)x Pi (cos 0) . Note that 
the polynomial Pi^"* (x, y) may now contain terms which are even in both x and y, yielding, in general, a non- vanishing 
contribution to the integral. To proceed, one thus has to be provided with the explicit form of Pq^''. 

The explicit form of the polynomial Pi^'' (x, y) is obtained by substituting for 5*0 and Si (and their gradients) from 
Eqs. ( pT| ) in Eq. ( p3b| ), taking 6r = St = and recalling Sx ^ x, Sy = y, and Uy = 0. One thereby obtains 

pW(:j,2;) = P^-)a;4 + P^-f)a;2y2+pfe)y4^ fora-t,r,x, (43a) 

P^^\x, y) = P^(-) x^y + P^^) xy\ (43b) 
where the various coefficients are explicitly given by 

P(-) = -\ [r'r'roi2ul - r^) + ro\2ut + 3uy, + 4)] , 

p('^)=-irg(l-rV2), (44) 
p(-) = -r,u,r{ul - rl/2\ P/^^^ = -rou,r{ul - rl), P/^^ = \rlu,r, (45) 
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Pi-) = 0, Pi"^) = \r,u,f{rl - 2ul), P^y^ = ir^u.f, 



(46) 



(47) 



In these expressions / = (1 — 2M/ro), r = u^, and all four-velocity components are evaluated at z. Note that the 
components P/^\ Pr^\ and Pi^^ consist of only terms which are even in both x and y. On the other hand, the y 
component P^^^"* contains only terms which are odd in both coordinates. 



Consider first the y component: Both terms oc x'^y and oc xy"^ of the polynomial P^^' yield, upon integrating, no 



KB 

contribution to Fy ' , and one immediately obtains 



P^(^) ^ 0. 



(48) 



l( 

The other components of Fa do not similarly vanish: Recalling x = pcosf and y = p sin if, and expressing eo in 

ay w 

^ro-^/»/^, (49) 
L PiicosO) 



the form eg = rQp{9){l + rg cos^ ip)^!"^ ^ we may write the double integral in Eq. ( ^8|) in the factorized form 



where 



27r 



■ dcosd, 



(50a) 



PA"'^ cos" (f + P^y^ cos^ LP sin> + PA^^ sin"* ip 



dip. 



(50b) 



(1 + Tq COS^ (p)5/2 

We now take advantage of the freedom we still have in specifying the function p{6), and make the convenient choice 

p = 2sin(6l/2). (51) 

With this choice, the integral becomes a standard one, reading simply 

h = (27r)-i (52) 

[see, e.g., Eq. (7.225-3) of The integral in turn, is a linear combination of standard elliptic integrals. It can 

be expressed as 



(53) 



where K[vj) and P(w) are two complete elliptic integrals of the first and second kinds, respectively, the argument w 
is given by 



rl + ul'' 



(54) 



and the coefficients Cq. and ha read 



with 



da = -iz{ro/u^fw-^'^ a(-)pl-) + af^^^P^"^) + a^^^P^^) 
ba = -^(ro/w,)«;-3/2 \b(-)p(-) + bi-y)p(-y) + 6to)pfe) 



a(^) ^ {w + 2){w-l), 



(55) 



-Av) - 



-2(w-l), 
(3w-2), 
= 2(w-l)2(w + l), 
6("^) = -(u;-2)(w-l), 
b^y^ = 2(1 - 2w). 



(56) 
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KB) 

The explicit form of the desired contribution Fa (for a = r,t, x) is finally obtained by inserting the values of 
Pa'K Pa'^\ and Pi^'' [given in Eqs. (|4^)-(|47|)] in the above expressions for and b^, constructing through Eq. 
@, and substituting in Eq. @. This yields 

UB) 1 utr[k{w) - 2E{w)] 

- 3 ^^37^ ' (57b) 



pKB) 1 r[K{^)-E{w)] 



where V =1 + u^/r^ 



KB) 

Note the remarkable fact that the contribution Fa is independent of I. 



D. Calculating 



Finally, let us evaluate f!^^\ Recalling F^"^^ — e^^Pi^^ and using Eqs. (23a) and (21a), Eq. ( ^9| ) becomes 

Fi^^) = -[L/(2^)](5°^ + UaUp)F'£, (58) 

where 

Pf= lim [ Sx'^e^^Pi{cos9)dn. (59) 

Note that we have already taken here the limit St 0, hence the integrand (cx Sx^) vanishes identically for j3 — t. 
Also, since cos9(p) and eo, given explicitly by 

eo = [rlp^ + + (urSr + u^x)^] , (60) 

are both even functions of y, the integral in Eq. (^9|) obviously vanishes for (3 = y. Hence, 

F^l = = 0. (61) 



Consider now Eq. ( |59| ) for the two remaining components, P — r,x. First, we change the integration variables to 
x,y. Since the Jacobian is d{6,(p)/d{x,y) = {pp')^^ (where p' = dp/d9), Eq. (|5^) becomes 



F^J' = lim / 5x^^€o^H{p)dxdy, (62) 

Si rO^ J 

where H{p) = Pi{cos9) sm9{pp')~^ . The function H{p) is a regular, even function of 9 (and of p), with H{Q) = 1. 
We thus write it as H{p) = 1 + p^H{p), where the function H{p) admits a regular (even) Taylor expansion at p = 0. 
Accordingly, we divide f!^ into two contributions, 

F!f = lun (I^+I^), (63) 

where 

/f = J Sx/^e^^dxdy, 

I2 = J Sx^€a^p^H{p)dxdy. (64) 
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Consider first the contribution : Near p = 0, the integrand in this term scales like Sr^ , thus the integrated singular 
contribution scales like Sr'^ . Hence, based on precisely the same argument applied in Appendix ^ with regard to the 
term Fjf"^ , we find that the integral is sufficiently regular to allow us interchanging the orders of the 6r ^ limit 
and integration: 

lim = I^{Sr = 0). (65) 
Doing so, we find that the contribution from to Pi^ vanishes for either jS = r or jS = x: For (3 = r, the integrand 



vanishes identically; for (3 = x, the integrand, evaluated at 5r = 0, becomes an odd function of x [see Eq. (41)], which 
vanishes upon integrating. 

To calculate the remaining contribution /f , we divide the domain of integration in Eq. ( |6^ into two regions: Let iJ'" 
denote the square —h < x,y < h, for some particular < h < 1 (say, h = 1/10); and let denote the remaining 

integration area over the sphere, outside i/"^. Correspondingly, we divide the integral /f into two contributions, as 
/f = /f + if Now, since the integrand of contains no singularity (the only singularity on the sphere occurs 
at a; = y = 0, which is located in if™), in evaluating lim^^^o* if°"'(<5?') we are allowed to interchange the limit and 
integration: 

lim /f™' =/f™V5r = 0). (66) 

Precisely as in the case of the integral considered above, this contribution is then found to vanish for either (3 = r 
or (3 = X. We are thus left with F!f = lim^r— o± -^f ™7 namely, 

= lim / / Sx^e^^dxdy. (67) 



-h J-h 



We proceed by considering separately the two components (3 = r and (3 — x. Let us begin with the r component: 
Re-scaling the integration variables as X = x / 5r and Y = y/5r^ we find 

/h/5r f-h j 8t 
\ \l±[X,Y)\-'^dXdY, (68) 
-hl&r J-h/6r 

where 

g± = e^/5r = ±[.9°, + rl{X'' + Y^) + («, + u^XfY'\ (69) 

and the ± sign refers to the sign of 5r. Note that €± (and hence the entire integrand) is independent of (5r, such that 
the 5r —> 0± limit becomes trivial: 

/oo poo 
/ [g^„ + rl{X^ + Y^) + {ur + u,Xf]-^/^dXdY = ±U. (70) 
-OO J — oo 

This is an elementary integral [see, e.g., Eq. (3.252-3) of [Q, in conjunction with Eq. (3.252-2) therein], yielding 

U = i2TT/ro)[rlul+glirl+ul)]-^/\ (71) 

Note the relation 

U = -^, (72) 

which stems directly from the "radial" geodesic equation of motion, (u' )^ — — {1 — u^/rDf. 

Consider next the case (3 = x. It is not possible to treat this case the same as the case (3 = r, hy changing 
the integration variables to X,Y: doing so, the integrand becomes Xe^'^, and the double- integral does not strictly 
converge at infinity. We therefore apply here a different method to evaluate the limit Sr 0^ in Eq. (^^. First, we 
express F^f as 



F'^ = lim / / xSa'^dxdy, (73) 

dr^O± J_h J-h 
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where, recall, So = = rf^ix^ + y^) + Qrr^r'^ + ("r'^'^ + Uxx)^ (with the limit t ^ to already taken). Now, So is 
quadratic in Sx" , and its derivative with respect to a; is a linear combination of both x and Sr. One easily obtains 
the relation 



where the coefficients a and (3 are given by 



X = aSo,x + PSr, 



f3- 



Substituting x from Eq. ( [74|) in Eq. j?^), we express as the sum of two integrals: 



plx 



lim 



a / 5*0,2; 5*0 dxdy + (3 / 5reQ dxdy 



ali + (3Iit 



(74) 



(75) 



(76) 



In what follows we show that 1^ vanishes, leaving us with only the contribution from la, which is just proportional 
to the r-component F!^ calculated above. 

Considering first li, we carry out the trivial integration over x, obtaining 



li = lim 



-1/2 



x—h 



dy. 



(77) 



The integration over y is then a standard one, but one do not need to carry it out explicitly: Observing that the 
integrand is now a regular function of y and 5r throughout the entire range of integration, we are allowed to interchange 
the limit and integration. Noticing then So{Sr = 0,x = +h) = So{Sr = 0, a; = —h), we immediately conclude 



/^ =0. 

Consider next la. Comparing with Eq. ( |67| ) (for the r components) we find simply la = F!^ , hence 

F^ = I3F![ ^ ±I3U. 

Having calculated all components of Pf, we may now construct F^^^^ through Eq. (psf). We obtain 



(78) 



(79) 



(80a) 



F'^^^ - T[L/{2tt)] [/-I + UriUr + f3Ux)] U^± 



(80b) 



f'^^^ = T[i/(2^)] [I3{rl + ul) + UxUr] U = 0, 



F. 



'(A) 

±y 



0, 



where we have substituted for U and /3 from Eqs. (|7^) and ([751), respectively. 

Note the remarkable fact that the contribution f!^^^ is precisely proportional to L. 



(80c) 
(80d) 



VIII. VALUES OF THE REGULARIZATION PARAMETERS 



In conclusion of the calculation carried out in the previous section, we have found that the /-mode direct force 
is composed of only two contributions: one — completely described by — is precisely proportional to L, and 

KB) 

the other — completely described by F^^ — is independent of L. No other powers of L are present. Recalling the 
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definition of the RP in Sec. |l^, we then conclude that the term -F^cf "* contributes only to the parameter and that 
the term i^^^"* contributes only to Ba- Recalling Eq. (|3^), we identify the RP as 

LA±^^q^F'^^\ B^ = q^F'^^\ C„=0. (81) 

Furthermore, from Eq. ( pj]) we immediately get Da = 0. The explicit values of A±a and Ba are then obtained by 

substituting the expressions derived above for the quantities F^^^'^'' — Eqs. (57), and (|^. 

To give a useful summary of the RP values thus obtained, we shall transform the angular coordinates a;, y back to 
the standard 6, (p coordinates, in which the orbit is equatorial (i.e., confined to 9 — tt/2). The quantities F^"{x) and 
F^"{x) are unaffected by this transformation, therefore the r and t components of all RP are unchanged. However, 
F^" and _F'^"' transform to Fg" and F^" in a manner which is not completely trivial, and we need to find the 
corresponding and (p components of the RP. Note that a-priori there is no guarantee that the RP will transform 
like vectors at the evaluation point, because the RP depend on F^" in the neighborhood of z; and the transformation 
(x, y) —^ (9, if) involves nontrivial functions of the angular coordinates, which may affect the mode decomposition. 
However, in Appendix ^ we show that in fact all the RP do transform (in this particular coordinate transformation) 
like four- vectors at z. It is trivial to show that at the evaluation point 

x.e = y.cp = x,^ = -yfi = 1 (82) 

[for concreteness we consider here the transformation described by a 7r/2 rotation about the horizontal axis = 0, 
which takes z from the pole to the point (6*0, <yJo) ~ (''■/2, —t^/'^) — see Fig. |l] and Appendix |^ for more details. Note, 
however, that the RP in the 0, ip coordinates do not depend on (/^g, due to the symmetry of rotations in if] Therefore, 
~ Bj. and Bg = —By (and the same for A^p, Ag). Note that u,j.{z) now becomes = C, the conserved azimuthal 
angular momentum. In the standard 9, ip coordinates, the RP are then given by Ag = Bg = 0, 



q"^ £ q^ r 



^±'- = Tl2— ' ^it^i^TT, ^^ = 0, (83a) 



'0 



where 



q^ £r[KH - mw)] . 
^* = 3 ^^757^ (83c) 



q^ r[KH - E{^)] 

Ca = Da^ 0, (83e) 
2, V=l + Cyrl (84) 



/ = (1 — 2M/ro), and = (u^)'^ — £^ — fV. Recall £ = —ut and C — u^p are the (conserved) specific energy and 
angular momentum parameters. 

We comment that the parameter A±a is normal to the four- velocity: A±au"' — 0. However, the parameter Ba (as 
the self force itself in our model) is, in general, not normal to u": An explicit calculation yields 

„ _ q^ rEjw) 

Finally, we give here the RP values for the special case of a radial geodesic, i.e., £ = 0: Noting, in this case, w — 0, 
V = 1, and -P = £^ — /, and recalling K{Q) = E{Q) = 7r/2, the nonvanishing components in Eqs. (p3|) reduce to 
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A^^f'^' = T^{£/f), ^S''"' = ±4 r, (86a) 



'^o '0 

^radial ^ ^ j-l (^2 _ ^ ^radial ^-f-gf. (86b) 

[The vanishing of iJ^^t^'a-i is obvious from symmetry considerations. Note that B^p vanishes at the hmit w despite 

the factor C in the denominator, because K{w) — E{w) = 0{w) = 0(£^).] These values are in agreement with the 
ones derived in [ p^ using the I mode Green's function expansion method. 

IX. CONCLUDING REMARKS 

The mode-sum scheme described by Eq. (|l2|), with the explicit RP values calculated in this paper, Eqs. (^3|), provide 
one with a practical means — yet one based on a physically well-established regularization scheme — for calculating the 
scalar self-force for any geodesic orbit around a Schwarzschild black hole. Recall that the full modes F^^"''^' needed 
for fully implementing this mode-sum scheme, are to be obtained from the /-modes of the scalar field which, in 
turn, are to be calculated using standard numerical techniques. 

The RP values derived here, Eqs. (^3|), were obtained independently by MNS using a different approach. In their 
analysis, MNS decomposed the direct field using the "standard" 9, ip coordinates (in which the motion is equatorial), 
in which case the rn 7^ modes contribute as well. MNS then derived an analytic expression for the contribution 
of each /,m mode of the direct force, expanded in powers of RI/7\ By explicitly summing up this expansion (and 
summing over m), MNS were able to recover all RP values. 

The RP values (p3) reduce, in the special cases of radial motion {C = 0) or circular motion (r = 0), to the values 
derived previously [l^] using a completely independent analytic approach (namely, by locally analyzing the /-mode 
Green's function, as we briefly mention in the introduction). In these cases, the values of the parameters Aa, Ba, and 
Ca have also been confirmed numerically, by calculating the full- force modes [ p^ . 

The RP calculation method presented in this paper is directly applicable to the more realistic case of the gravitational 
self-force acting on a mass particle, as well as to the case of the electromagnetic self-force acting on an electrically- 
charged particle. Both cases shall be treated in an accompanying paper where we obtain the gravitational and 
electromagnetic RP for general orbits in Schwarzschild spacetime (the results in the gravitational case were provided 
in [^). The extension of our scalar field analysis to the gravitational and electromagnetic cases involves several 
complexities which require a special care. In particular, one has to tackle the technical issue of extending the four- 
velocity vector off the worldline |0] . A more fundamental issue concerns the gauge dependence of the gravitational 
self- force 01. 
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APPENDIX A: DERIVATION OF So AND Si 

In this appendix we calculate the two leading terms in the expansion of S* = S{x, z) [the square of the geodesic 
distance from the point x to the geodesic z{t)] in powers of 5x^ = x^ — . This expansion takes the form 

5 = 5o + 5i + 52 + • • • , (Al) 

in which the term Sn is of homogeneous order 5a;"'+^, and we wish to calculate So and Si. 

In flat space, using Cartesian coordinates y" (with j/" = at z), we obviously have S = {rjap + UaUi3)y°'y^ = Sq, 
where ria/3 is the flat space metric. In curved space (or in curvilinear coordinates), each of the terms S„ (like S itself) 
will be a certain function of gafj and its derivatives. From simple dimensionality considerations, it is clear that Sq 
may not include any derivatives of ga/3, and that 6*1 may include only first-order derivatives of the latter (in addition 
to gap itself) . 
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Let y" be locally-Cartesian coordinates at the evaluation point z, with = at z. Namely, at a; = z, the metric 
functions in the coordinates are just t]^^, and their first-order derivatives vanish. Since no second or higher-order 
derivatives appear in S up to the desired order, we must have 



S=ir^^p+u'y0)y'^y^ + Oiy'), 



(A2) 



where a prime denotes vectorial components in the y" coordinate system. We now transform from back to our 
original coordinates x^^. Recall that S* is a bi-scalar, and is hence invariant under this transformation. Writing the 
Taylor expansion of y"((5x^). 



dy"' 1 d^y" 



2 dxt^dx" 



Sx^'Sx" + 0{5x^) 



(A3) 



(in which all coefficients are evaluated at z), and substituting it in the right-hand side of Eq. (A2), we find 



(^./3 + u„u^)^^ 



5x^5x^5x'' +0{5x'^). 



Comparing this to Eq. ( Al), we identify the first and second terms on the right-hand side with and 5*1, respectively. 
Using the obvious tensorial transformation rule, we find 



■5*0 = {gxe + uxu^)5x^5x^. 
To calculate we need the second-order transformation coefficients, which are given by 



(A4) 



^x^'^x^' dx 



[see, e.g., Eq. (3.2.11) of [||]. Therefore, 

^1 - 



, , dy°' dy^ 

(gXe + UxUe)^ Sx^Sx'^dx". 



Sx^Sx^'Sx" 



Recalling that 



we finally obtain 



gxeTlJx^Sx^Sx" = -g^.^xSx^Sxf'Sx-', 
Si = {uxu^T^1,+gl0^^/2) Sx'^Sx^Sx^ 



(A5) 



APPENDIX B: INTERCHANGEABILITY OF THE r ^ n, LIMIT AND THE LEGENDRE INTEGRAL 

In this appendix we explore the interchangeability of the limit and integration in Eq. (|36| ) — an issue crucial for the 
calculation carried out in Sec. VII. For convenience, let us write F^^^ = [L / {2tt)]F^^\ where 



F. 



{W)l 
±a 



lim 



2tt 

d9 I d(pF^^\r, to, 0,(p)Pi{cos 9) sine. 





Here, recall, W stands for A, B, or C, with 

-3p(l) p(S)^,-5p(4) p(C)^ -7p(7) 



(Bl) 



(B2) 
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where eo is given explicitly in Eq. (^^, and Pi"-* represents a polyn omi al of homogeneous order n in 5x'^ = — z^. 
We shall show that interchanging the limit and integration in Eq. (Bl) is valid for W — B,C, and explain why our 
proof fails in the case W = A. 

We begin by considering the case W = B. Let R= \Sr\^ for some < s < 1. Consider the range of small Sr (such 
that R < 1). We spht the ^-integral in Eq. (Bl) into three domains: (i) 9 < R, (ii) R < 9 < 1, and (iii) 1 < 9 < n. 
Correspondingly, the above double-integral can be expressed as 2i + Tu + Tm , hence 

lim I^{5r)+ lim Iu{5r) + lim I^^^{5r). (B3) 

5r^0± f)>^0± <5r^0± 

Consider first the internal integral Xi : 

I,{Sr)^ d9 dipFi^'>{5r,9,(p)P\cos9)sm9. (B4) 

^0 "'0 

Since at the relevant limit i? — + 0, we have 9^1 throughout the range of integration, and we may approximate this 
integral [using siii9d9 ~ pdp, as well as P\cos9) ~ 1 and p{R) ~ R] as 

2", {Sr) ~ r pdp r d^ Fi^) ((5r, p, ^) . (B5) 
Jo Jo 

Now, Fi^^ = eQ^Pi'^K We may bound eo and Pi^'' as eo > ci\Sr\ and |-Pi^^| < C2/9^, where hereafter c„ are some 
positive constants. Consequently, Fa^^ can be bounded as < C3\Sr\^^p'^ < C3\Sr\~^R'^. Since the "integration 

area" in Eq. ( p35| ) is 7r_R^, we then obtain the upper bound 

< nc3\Sr\-^R'^ = Ci\Srf''^ . (B6) 

Taking, e.g., s = 0.9, we find that \Xi\ < C4|5r|" '* and hence it vanishes at the limit 6r 0^: 

lim J,{dr) = 0. (B7) 

dr-^Q± 

Consider next Ta , which is defined by the double- integral 

I,,{Sr)= d9 dipFi^\5r,9,ip)P\cos9)sin9. (B8) 

JR. Jo 

( B) ( ( B) 

Let us introduce the quantity AFa = Fa {6r) — Fa [Sr — 0), and the corresponding integral 

AI,,{6r)= d9 d(pAFj^^'>{dr,9,(p)P\cos9)sm9, (B9) 

JR. Jo 

such that 

r2TT 



1ii{5r) 



pi p2TT 

AIi,{6r)+ d9 dipF^^^Sr ^0,9, 'fi)P' {cos 9) sin 9. (BIO) 

JR JO 



Since s < 1, in the entire range (ii) we have \Sr\ <ti R < 9 (for \Sr\ <C 1), and thus |^7'| <C p. We may approximate 
AFi^^ at |(5?i < 1 as 

AF^) ^ d(eo ^p(4))/d(fr) X Sr. (Bll) 

dr—O 



From Eq. ( |60| ) (and recalling \x\ < p) we observe that, at 6r = 0, ep ^ and d{eQ^)/d{6r) are bounded from above 

by cx p^^ and oc const, respectively. At Sr = we may also upper bound P^^ and dP^^ /d{Sr) by oc and (x p^, 
respectively. Applying these bounds in Eq. ( |B11| ) we then obtain 

\AFi^^\<c,\Sr\p-'. (B12) 

[Note that in Eq. ( |B11| ) we have assumed that either S or P^"^ (or both) include terms linear in Sr. In special cases 

( R) 

where only terms quadratic in Sr are present in both quantities, we shall instead arrive at the bound \AFa \ < 
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ceSr'^p Since \Sr\/p <ti 1, this is even smaller than the bound ( B12| ), and so the entire derivation below remains 
vahd.] 

Using Eq. (pl|) and recalhng |P'(cos6')| < 1, we can now bound AIii{5r) by 



\/Slu{5r)\ < C5 I del dip sin e \6r\p-^ =2ttc5 dOsine 
Jr Jo Jr 



\5r\p' 



Since siii9/p and dO/dp are bounded in domain (ii), the last integral can now be bounded as 

fp(i) 



\AIu{Sr)\ < ce\5r\ f p-'dp. 

Jp(R) 



(B13) 



(B14) 



[We point out here that in the way we have chosen p{6) in Sec. [VI I| — see Eq. ( ^l| ) — dd/dp is unbounded &t 6 ^ tt. It 
is this divergence that forced us to terminate domain (ii) at 6* = 1, and to introduce domain (iii).] Upon integration 
we find 



|AI,,(fr)| <C6|fr|[lnp(l)-lnp(i?)]. 
For small R we have p{R) ~ R— |(5r|^; therefore 

\AIuiSr)\ < C6lnp(l)|fr| -sc6|<5r|ln|<5r|. 
Clearly, both terms vanish as 5r — > 0^, hence 

lim AJu(5r) = 0. 

<5r^0± 



(B15) 
(B16) 
(B17) 



From Eq. (BIC) we then have 

lim Ir, {5r) = lim / dO dip F^^^ {6r = Q,0,p)P^ (cos 9) sin 9 



pi p2-K 

d9 dipF^^^Sr = 0,9, Lp)P\cos 9) sin 9. 
Jo Jq 



Finally, consider the third contribution, Im, defined by 



I,u{Sr)= d9 dp Fi^\Sr, 9, ip)P\cos 9) sin 9. 



In full analogy with the analysis above, we define 



AI^u{Sr)= d9 dpAFi^\5r,9,p)P\cos9)sin9. 
J I Jo 



(B18) 



(B19) 



(B20) 



The above bound, \AFi^^\ < C5\Sr\p-^, is valid in range (iii) too. Since in this range p is bounded from below, we 

may now write \AFa \ < C7\Sr\; and since P {cos9) and sin 9 arc both bounded by unity, the entire integrand is then 
bounded by C7|<5r|: 



d9 J dp = 27T^C7\6r 



(B21) 



Again, this quantit y van ishes at the limit Sr 0; hence, by the same considerations used above for range (ii) [see the 
chain of equations (|B18| )], we obtain 



p7T p'2'K 

lim I^.u{Sr)= d9 dpFi^\Sr = 0,9, p) P^ {cos Vjsmt 



(B22) 



Substituting Eqs. (B7), (B18), and (B22) in Eq. (B3), we obtain 
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/•TT p'2i'K 

F^^^^ = \ de dipFi^^Sr ^0,9, f)P\cos 6) sine. (B23) 
Jo Jo 

Namely, in the calculation of F_^^ ' — and thus also F^^ ' — we are allowed to interchange the limit 5r ^ and the 

integration. Note that since Fa admits a well-defined limit at Sr ^ (except at 6* = — which, however, was shown 
not to affect the integral), we have omitted its ± label. 

The same proof can immediately be applied to f'^^\ Evaluating Ji, we find this time that IFq'^'I is bounded from 
above by cslSrl'"^ < cs\5r\~'^ R'^ , hence (taking again s = 0.9) 

< cs\Sr\-'R^ = cs\Srf'-' = cg\Sr\'-^ ^ (B24) 
as dr 0. Evaluating next Ala, we obtain this time [in analogy with Eq. ( B12 )] lAFi*^-*! < cio\Sr\/ p, and hence 



/•p(i) 

AJ,,(fr)| < ciil^rl / p"dp<ci2\Sr\, (B25) 

Jp{R) 



which again vanishes as (5r ^ 0. The calculation of AJm proceeds exactly as for = i? (the only difference is that 

now AF^^'' (X \Sr\/p instead of |(5r|//9^, but this does not affect the above evaluation of AZm in any way). Again we 
find |AXMi((5r)| < 27r^ci3|(5r|, which vanishes at the limit Sr —> 0. We conclude that the limit and integration may be 
interchanged for = C as well: 

F^^^^ = dd dip fP {6r ^Q,0,ip)P' (cos d) sin 9. (B26) 
Jo Jo 

Finally, it is instructive to see how the above type of arguments fail for W = A. Since F^^^ = eQ^Pi^\ in evaluating 
Xi one obtains \Xi\ < ci4\Sr\~^R^ = cii\Sr\'^^^~^\ Then, evaluating AJu, one obtains 

\AIu{5r)\ <ci5\Sr\ / p'^dp, (B27) 

Jp{R) 

which at the limit of small R yields \AIii{Sr)\ < ci5\Sr\/R = cisj^rj^"*. Obviously, for any s < 1 the bound for Ti 
will fail to vanish as dr 0, and for any s > 1 the bound for Ala will fail to vanish at this limit. [Note also that in 
the case s > 1 the inequality \5r\ <C p, used above in evaluating Ala, is no longer valid throughout range (ii)]. In fact, 
it becomes evident from the explicit calculation in Sec. VII, that for W — A the limit Sr Q cannot be interchanged 
with the integration. 



APPENDIX C: TRANSFORMING TO EQUATORIAL ORBIT 



In Sec. VII we analyzed the multipole decomposition of the direct force in a Schwarzschild coordinate system in 
which the particle is momentarily at the pole. We then transformed to locally Cartesian angular coordinates x, y and 
calculated the RP in the system = {t,r,x,y). Usually (e.g., in numerical calculations) one adopts a more natural 
pair of angular coordinates, in which the particle's orbit is confined to the equatorial plain (9 = 7r/2) — throughout 
this appendix we shall denote this system by x^^ = (t, r, 9, (p). The goal of this Appendix is, given the RP values in 
the system xf^^, to obtain the corresponding values in the system a;^q. 

In the system x'^^ we have = and = £, where, recall, C is the conserved specific angular momentum. Let the 
particle be momentarily at {9,ip) — {9o,po) ~ ('''/2, — 7r/2).[| Consider the set of spherical coordinates 9,ip, defined 
such that the particle is located at their pole, 9 = 0, and ip = Q,7r coincides with 9 = tt/2 — see Fig. 0. (These are, 
in fact, the same spherical coordinates used throughout the paper; here we merely use a different notation, as the 



^ We consider liere a specific value of (po in order to simplify the following expressions, and to make the correspondence 
between the x,y and x,y coordinates (see below) easily apparent. Note, however, that our final result — the RP values in the 
system x^fq — do not depend on the choice of ipo, as the Z-mode decomposition of the direct force in that system is invariant 
under rotations about the polar axis. 
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symbols 0, (p are reserved for the angular coordinates of the x^^ system.) The "locally Cartesian" coordinates at z are 
then given by [see Eq. (31)] 

X = p{9) cos = 2 sin(0/2) cos (p, 

?/ = p(^)sin(^ = 2sin(0/2)sin(^. (CI) 

Relating the spherical coordinates 0, (p to the standard pair O^Lp is a, straightforward geometrical problem, and one 
finds 

cos^ = — sin6'sin(^, cot = tan 6* cos (p. (C2) 
This allows us to express x, y directly as functions of 0, Lp. We obtain 

Vl — sm6'sin(^ 

y{e,^)^2'/^ '""'l . , (C3) 

V 1 — sm y sm ip 

which describes explicitly the transformation between and x'^^. Note that this transformation is regular on the 
entire sphere [except at the point {9-,^) = (7r/2,7r/2), which, however, is irrelevant for our analysis]. 




FIG. 1. A Sketch showing the various quantities involved in constructing the coordinate transformation (x,y) — > {d,ip). Shown 
is the "northern" hemisphere t = r = const, 6 > 0. The particle moves along an equatorial orbit, and is momentarily located at 
(6, f) = (7r/2, — 7r/2). x,y is a pair of "locally Cartesian angular coordinates" at the particle's location, as described in the text. 

Given the above transformation rule, the desired components Fg" and F^f"' are constructed as 



tg - Xfit^ 
pdir _ pdir 



(C4) 



It is now useful to consider the Taylor expansion of the various partial derivatives about the particle's location: 
Introducing A6 = — 7r/2 and Ap = ip — (— 7r/2), we find 



x.e 



-AdAip + 0{Sx^), 



8 ^ 8 



0{Sx^), 



y,e = -l-lA^' + lA9^ + OiSx'), 



2/,^ - -]A0Aip + OiSx^), 



(C5a) 
(C5b) 
(C5c) 
(C5d) 
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where O(fe^) represents corrections of fourth order in A6' and A(^. 
obtain, near the particle's location, 



Substituting these exapnsions in Eqs. (C4), we 



(C6) 



where corrections are due to terms of the form eg Pg "' (re call the notation introduced in Sec. ^ ), and higher-order 
terms that vanish at x z. From the analysis of Sec. VII it is clear that such correction terms do not contribute to 
neither i^i"^'' nor F^\ and their contribution to F^^ vanishes in the multipole decomposition. Hen ce, none of the 



RP will be affected by omitting these correction terms and replacing the approximation in Eq. (C6) with an exact 



equality. Consequently, we find that the RP transform under x^^ 



as vectors at z, namely 



Re 



Rip — Rx 5 



(C7) 



where R^ stands for any of the RP (obviously, Rt and Rr do not change). Note also the relations ug = —Uy{z){= 0) 
and ~ Ux{z) {ut and do not change). 



[1] Up to date infomiatat ion about the joint NASA-ESA project LISA can be found in the project's webpage at 



http:/ /lisa.jpl.nasa.gov/ 
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